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Abstract: In this paper we consider a continuous-time anisotropic swarm model 
with an attraction/repulsion function and study its aggregation properties. It is 
shown that the swarm members will aggregate and eventually form a cohesive cluster 
of finite size around the swarm center. We also study the swarm cohesiveness when 
the motion of each agent is a combination of the inter-individual interactions and 
the interaction of the agent with external environment. Moreover, we extend our 
results to more general attraction/repulsion functions. The model in this paper is 
more general than isotropic swarms and our results provide further insight into the 
effect of the interaction pattern on individual motion in a swarm system. 



u 

6 
cr 



> 

oo 



Keywords: Autonomous mobile agents, biological systems, multi-agent systems, 

VO ' swarm intelligence, aggregation. 

O ' 

1 Introduction 
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In nature swarming can be found in many organisms ranging from simple bacteria to more 
advanced mammals. Examples of swarms include flocks of birds, schools of fish, herds 
of animals, and colonies of bacteria. Such collective behavior has certain advantages 
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rj such as avoiding predators and increasing the chance of finding food. Recently, there has 

been a growing interest in biomimicry of forging and swarming for using in engineering 
applications such as optimization, robotics, military applications and autonomous air 
vehicle PQ [I|. Modeling and exploring the collective dynamics has become an important 
issue and many investigations have been carried out [H] JH]- However, results on the 
anisotropic swarms are relatively few. The study of anisotropic swarms is very difficult 
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though anisotropic swarming is a ubiquitous phenomenon, including natural phenomena 
and social phenomena. 

Gazi and Passino [2] proposed an isotropic swarm model and studied its aggregation, 
cohesion and stability properties. Subsequently, Chu, Wang and Chen jlj generalized 
their model, considering an anisotropic swarm model, and obtained the properties of 
aggregation, cohesion and completely stability. The coupling matrix considered in [4J is 
symmetric, that is, the interaction between two individuals is reciprocal. In this paper, 
we will study the behavior of anisotropic swarms when the coupling matrix is completely 
asymmetric. The results given in this paper extend the corresponding results on isotropic 
swarms [2j and anisotropic swarms [I] to more general cases and further illustrate the effect 
of the interaction pattern on individual motion in swarm systems. Moreover, we also study 
the aggregation properties of the anisotropic swarm under an attractant /repellent profile. 

In the next section we specify an "individual-based" continuous-time anisotropic swarm 
model in Euclidean space which includes the isotropic model of [2] as a special case, and we 
also study the agent motion when the external attractant/repellent profile is considered. 
In Section 3, under some assumption on the coupling matrix, we show that the swarm 
exhibits aggregation. In Section 4, we extend the results in Section 3 by considering a 
more general attraction/repulsion function. We summarize our results in Section 5. 

2 Anisotropic Swarms 

We consider a swarm of iV individuals (members) in an n-dimensional Euclidean space. 
We model the individuals as points and ignore their dimensions. The equation of motion 
of individual % is given by 

N 

i^X^/Ca:*-^), z = l,...,JV, (1) 

where x l G R n represents the position of individual i; W — [wij] € R NxN with Wij > for 
all i,j = 1, • ■ ■ ,N is the coupling matrix; /(■) represents the function of attraction and 
repulsion between the members. In other words, the direction and magnitude of motion 
of each member is determined as a weighted sum of the attraction and repulsion of all the 
other members on this member. The attraction/repulsion function that we consider is 

f(y) = -y(a-bex P (-^Pj\ (2) 

where a, b, and c are positive constants such that b > a and ||y|| is the Euclidean norm 
given by \\y\\ = y/y^y. 

In the discussion to follow, we always assume wa = 0,i = 1, • • • ,N in model (JIJ. 
Moreover, we assume that there are no isolated clusters in the swarm, that is, W + W T 
is irreducible. 



Note that the function /(•) is the social potential function that governs the interindi- 
vidual interactions and is attractive for large distances and repulsive for small distances. 
By equating f(y) = 0, one can find that /(•) switches sign at the set of points defined 
as J = [y = Oor ||y|| = 5 = yc In (b/a) } . The distance 5 is the distance at which 
the attraction and repulsion balance. Such a distance in biological swarms exists indeed 
j3]. Note that it is natural as well as reasonable to require that any two different swarm 
members could not occupy the same position at the same time. 

Remark 1: The anisotropic swarm model given here includes the isotropic model of 
|2] as a special case. Obviously, the present model ((H) can better reflect the asymmetry 
of social, economic and psychological phenomena [Tlj-|19|. 

In the above model, the agent motion was driven solely by the interaction pattern 
between the swarm members, i.e., we didn't consider the external environment's effect on 
agent motion. In what follows, we will consider the external attractant/repellent profile 
and propose a new model. 

Following ^T|, we consider the attractant/repellent profile a : R n — > R, which can 
be a profile of nutrients or some attractant/repellent substances (e.g. nutrients or toxic 
chemicals). We also assume that the areas that are minimum points are "favorable" to the 
individuals in the swarm. For example, we can assume that a(y) > represents a noxious 
environment, a(y) = represents a neutral, and cr(y) < represents attractant or nutrient 
rich environment at y. (Note that a(-) can be a combination of several attractant/repellent 
profiles) . 

In the new model, the equation of motion for individual i is given by 

N 

x 1 = -hiW x ra(x l ) + ^2 w ijf( xl ~ xJ )> i = 1> " ' ' >-W> ( 3 ) 

3=1 

where the attraction/repulsion function /(■) is same as given in (J2J), hi G R + = (0, oo), 
and Wij is defined as before. — hiW x i(j{x % ) represents the motion of the individuals toward 
regions with higher nutrient concentration and away from regions with high concentration 
of toxic substances. We assume that the individuals know the gradient of the profile at 
their positions. 

In the discussion to follow, we will need the the concept of weight balance condition 
defined below: 

Weight Balance Condition : consider the coupling matrix W = [w^] G R NxN , for 
all i = 1, • • • , N, we assume that ^, =1 w%j = J2j=i w ji- 

The weight balance condition has a graphical interpretation: consider the directed 
graph associated with the coupling matrix, weight balance means that, for any node in 
this graph, the weight sum of all incoming edges equals the weight sum of all outgoing 
edges |5 a . The weight balance condition can find physical interpretations in engineering 
systems such as water flow, electrical current, and traffic systems. 



3 Swarm Aggregation 



In this section, theoretic results concerning aggregation and cohesiveness of the swarms 
(JTJ) and (J3J) will be presented. First, it is of interest to investigate collective behavior of 
the entire system rather than to ascertain detailed behavior of each individual. Second, 
due to complex interactions among the multi-agents, it is usually very difficult or even 
impossible to study the specific behavior of each agent. 

~k Sj=i x% i an d denote fy = exp ( — 



Define the center of the swarm members as x 
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) . We first consider the swarm in Q, then the equation of motion of the center is 
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If the coupling matrix W is symmetric, by the symmetry of /(•) with respect to the origin, 
the center x will be stationary for all t, and the swarm described by Eqs. (P) and Q will 
not be drifting on average [1]. Note, however, that the swarm members may still have 
relative motions with respect to the center while the center itself stays stationary. On the 
other hand, if the coupling matrix W is asymmetric, the center x may not be stationary. 
An interesting issue is whether the members will form a cohesive cluster and which point 
they will move around. We will deal with this issue in the following theorem. 

Theorem 1: Consider the swarm in (^) with an attraction/replusion function /(•) in 
(0). Under the weight balance condition, all agents will eventually enter into and remain 
in the bounded region 
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and M = Yl w if ^2 denotes the second smallest real eigenvalue of the matrix L + L 
with 
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Q provides a bound on the maximum ultimate swarm size. 
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Proof. Let e l = x l — x. By the definition of the center x of the swarm and the weight 
balance condition, we have 
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Then, we have 
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To estimate e\ let V = Y2i=i ^ ^e the Lyapunov function for the swarm, where 
Vi = \e lT e % . Evaluating its time derivative along the solution of system (jljl. we have 
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where e = (e , ■ • • , e ) , L £g> I is the Kronecker product of L and 7 with L as defined 
in Eq. (J3J) and 7 the identity matrix of order n. 

Note that each of the functions exp ( — " x ~ x " ) ||x* — x- 7 !) is a bounded function whose 
maximum is achieved at || x 1 —x° \\ = v/c/2 and is given by -y/c/2exp(— (1/2)). Substituting 
this into the above inequality and using the fact that ||e J || < V2V, we obtain 

V < -ae T (L ® J)e + 2bMy/cexp ( - -) V* . (5) 

To get further estimate of V', we only need to estimate the term e T (L ® I)e. Since 

e T {L ® J)e = ^e T ((L + L T ) ® j)e, 



we need to analyze e T ((L + -7 T ) ® I)e. First, consider the matrix L + L T with L defined 



in Eq. (jlj), we have L + L T = [l^], where 
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(6) 



Under the weight balance condition, we can easily see that A = is an eigenvalue of 
L + L T and u = (I,- ■ ■ , l) T with / ^ is the associated eigenvector. Moreover, since L + L T 
is symmetric and W + W T (hence, L + L T ) is irreducible, it follows from matrix theory JSJ 
that A = is a simple eigenvalue and all the rest eigenvalues of L+L T are real and positive. 
Therefore, we can order the eigenvalues of L + L T as = Ai < A 2 < A3 < ■ ■ ■ < A n . 



Moreover, it is known that the identity matrix / has an eigenvalue \x 
and n linearly independent eigenvectors 
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By matrix theory j5j, the eigenvalues of (L + L T ) <g) I are Aj/i = A« (of n multiplicity 
for each i). Next, we consider the matrix (L + L T ) <g) I. A = is an eigenvalue of n 
multiplicity and the associated eigenvectors are 
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Therefore , e T ((L + L T ) <g> /)e = implies that e must lie in the eigenspace of (L + L T ) ®I 
spanned by eigenvectors v l , ■ ■ ■ ,v n corresponding to the zero eigenvalue, that is, e 1 = 
e 2 = ■ • • = e N . This occurs only when e 1 = e 2 = ■ • ■ = e N = 0. However, this is 
impossible to happen for the swarm system under consideration, because it implies that 
the N individuals occupy the same position at the same time. Hence, for any solution 
x of system (JTJ), e must be in the subspace spanned by eigenvectors of (L + L T ) <g> / 
corresponding to the nonzero eigenvalues. Hence, e T ((L + L T ) ® l)e > A 2 ||e|| 2 = 2\ 2 V. 
From ©, we have 

V < -a\ 2 V + 26Mv^exp(-|)\/^ 

aX 2 V 1 / 2 -2bM y /dexp(- 
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Therefore, any solution of system (PJ will eventually enter into and remain in Q. 
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Remark 2: The discussions above explicitly show the effect of the coupling matrix 
W on aggregation and cohesion of the swarm. 

Remark 3: The weight balance condition is more general than the case when the 
coupling matrix W is a symmetric matrix E] • 

Remark 4: Theorem 1 shows that the swarm members will aggregate and form a 
bounded cluster around the swarm center. 

Remark 5: From Theorem 1, we see that, under the weight balance condition, the 
motion of the swarm center only depends on the repulsion between the swarm members. 

From the above discussions, we know that if we ignore the influence on agent mo- 
tion from external environment, under the weight balance condition, the motion of the 
swarm center only depends on the repulsion between the swarm members, and all agents 
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will eventually enter into and remain in a bounded region around the swarm center. In 
what follows, we will study the aggregation properties of the swarm system when the 
attractant/repellent profile is taken into account. 

The equation of the motion of the swarm center now becomes 



N , N 
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Before we discuss cohesiveness of the swarm, we first make an assumption. 
Assumption 1: There exists a constant a > such that 

II V y a(y) || < o 7 , for ally. 

Assumption 1 implies that the gradient of the profile is bounded. This assumption is 
reasonable since almost all profiles we encounter such as plane and Gaussian profiles are 
with bounded gradient. 

The following theorem shows that the swarm system still exhibits aggregation behavior 
when the external profile is taken into account. 

Theorem 2: Consider the swarm in Q with an attraction/replusion function /(•) 
in @. Under the weight balance condition and Assumption 1, all agents will eventually 
enter into and remain in the bounded region 

N 
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where 
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and M and A2 are defined as in Theorem 1. Q provides a bound on the maximum ultimate 
swarm size. 

Proof. Let e l = x % — x. By the definition of the center x of the swarm and the weight 
balance condition, we have 
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Define the Lyapunov function as V = Yli=i Vh where Vi = |e e\ Evaluating its time 



derivative along solution of the system ©, we have 
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Furthermore, by assumption, we have 
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By analogous discussions as in the proof of Theorem 1 , we have 
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Therefore, any solution of system @ will eventually enter into and remain in Q. 
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Remark 6: Theorem 2 shows that, with bounded attractant/repellent profile, the 
swarm members will aggregate and form a bounded cluster around the swarm center. 
The motion of the swarm center depends on the repulsion between the swarm members 
and the weighted average of the gradient of the profile evaluated at the current positions 
of the individuals. 

Of course, not all the profiles are bounded. In the case of unbounded profile, in order 
to ensure the swarm to be ultimately bounded, the gradient of the profile at x % should 
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have a " sufficiently large" component along e l so that the influence of the profile does not 
affect swarm cohesion. The following theorem addresses this issue. 

Theorem 3: Consider the swarm in (J3J) with an attraction/replusion function /(•) in 
(J2J). Assume that there exist constants A l a , % = 1, • • • ,N, with A a = mmA l a > — ^ such 

i " 

that 
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for all x l and x k . Then, under the weight balance condition, all agents will eventually 
enter into and remain in the bounded region 
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i=l 



where 
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2bMV2cexp(-±] 



a\ 2 + 2A a 

and M and A2 are defined as in Theorem 1. Q provides a bound on the maximum ultimate 
swarm size. 

Proof. Following the proof of Theorem 2, from (J7J), we have 
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By analogous discussions as in the proof of Theorem 1 , we have 

V < -(a\ 2 + 2A a )V + 26M v ^exp(-|)y 1 / 2 
{a\ 2 + 2A a )V 1 / 2 - 26Mv^exp(-i 
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Therefore, any solution of system (J3J) will eventually enter into and remain in Q. D 



4 Further Extensions 

In Sections 2 and 3 we considered a specific attraction/repulsion function f(y) as defined 
in @. In this section, we will consider a more general attraction/repulsion function f(y). 
Here f(y) is still the social potential function that governs the interindividual interactions 
and is assumed to have a long range attraction and short range repulsion nature. Following 
[10 , we make the following assumptions on the social potential function: 
Assumption 2. The attraction/repulsion function /(•) is of the form 

f(y) = -y[fa(\\y\\)-fr(\\y\\)UeR\ (8) 

where f a : R + — > R + represents (the magnitude of) attraction term and has a long 
range, whereas f r : R + — > R + represents (the magnitude of) repulsion term and has a 
short range, and R + stands for the set of nonnegative real numbers, \\y\\ = \/y T y is the 
Euclidean norm. 

Assumption 3. There are positive constants a, b such that for any y e R n , 

fa(\\y\\)=a } f r {\\y\\) < ±r. (9) 

\\y\\ 

That is, we assume a fixed linear attraction function and a bounded repulsion function. 

Analogous to Theorems 1-3, in this case, we can also obtain the following three theo- 
rems. 

Theorem 4: Consider the swarm in ((TJ) with an attraction/replusion function /(•) in 
(jHJ) and 0. Under the weight balance condition, all agents will eventually enter into and 
remain in the bounded region 

r N 



< x : ^T) ||x l - x\\ 2 < p 2 >, 

^ i=i ' 



where p = -^— ; and A2 and M are defined as in Theorem 1; Vt* provides a bound on the 
maximum ultimate swarm size. 

Theorem 5: Consider the swarm in Q with an attraction/replusion function /(•) 
in (jHj) and 0. Under the weight balance condition and Assumption 1, all agents will 
eventually enter into and remain in the bounded region 

N 



i=l 

where 



^ = {*:|>*-x|| 2 <p 2 }, 
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_ 4jM + 4o(EtA) . 

and M and A2 are defined as in Theorem 1. Q provides a bound on the maximum 
ultimate swarm size. 
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Theorem 6: Consider the swarm in (JHJ) with an attraction/replusion function /(•) in 
and (jUJ)- Assume that there exist constants A l a , i = 1, • • • ,N, with A a = minA^ > 



■sf- such that 
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e iT 



> A 1 IIpMI 2 



hiV^a^) - — Y^ h k V x ka(x k ) 
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for all x % and x k . Then, under the weight balance condition, all agents will eventually 

enter into and remain in the bounded region 

N 



»=1 

where 



H* = jx:^||x 4 -x|| 2 <A 



AbM 
P 



a\ 2 + 2A a ' 

and M and A2 are defined as in Theorem 1. Q provides a bound on the maximum 
ultimate swarm size. 

Following the proof of Theorems 1-3, we can prove Theorems 4-6 analogously. 

5 Conclusions 

In this paper, we have considered an anisotropic swarm model and analyzed its aggrega- 
tion. Under the weight balance condition, we show that the swarm members will aggregate 
and eventually form a cohesive cluster of finite size around the swarm center. The model 
given here is a generalization of the models in |2], jl], and jTTj, and can better reflect the 
asymmetry of social, economic and psychological phenomena [T4*]-|19j. 
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